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1  Summary 


We  carried  out  numerical  simulations  on  multilayer  plasmon  excitations  by  resonant  energy 
transfer.  A  formal  derivation  of  the  transfer  rate  from  the  plasmons  to  the  Wannier  excitons  was 
performed.  We  stimulated  fluorescence  and  compared  our  results  with  the  conventional 
organic/inorganic  hybrid  organic  light-emitting  diode  (OLEDs). 

2  Introduction 

We  formulated  electron  energy  loss  spectroscopy  (EELS)  for  multilayer  epitaxial  grapheme. 

Our  derivation  of  the  general  angles-resolved  EELS  spectra  for  transmission  and  reflection  is 
based  on  the  prescribed  electron  trajectory.  In  particular,  we  were  interested  in  application  of  the 
theory  to  various  multi-layer  electron  gas  (MLEG)  configurations,  i.e.,  epitaxial  vs.  exfoliated 
graphene. 

We  did  a  comparison  with  the  available  experimental  data  to  see  the  effect  of  the  gap  between 
the  valence  and  conduction  bands. 

3  Methods,  Assumptions,  and  Procedures 

We  initiated  a  study  of  the  plasma  instability  in  graphene  double-layers  in  interaction  with  a 
thick  conductor  when  an  electric  current  is  passed  through  one  of  the  layers. 

For  optoelectronic  detectors,  we  investigated  the  relation  between  the  induced  and  driving 
currents  via  nonlocal  response. 

Comparison  of  the  plasmon  generation  with  conventional  two-dimensional  electron  gas  (2DEG) 
configurations  was  initiated.  The  role  of  the  gap  and  interband  polarization  was  investigated.  The 
aim  was  converting  plasmonic  energy  to  THz  radiation  using  a  grating. 

We  carried  out  a  study  of  dressed  Dirac  electrons  trapped  in  a  quantum  dot  by  coherent  nonlinear 
THz  spectroscopy. 

In  our  published  paper  “Tunable  electromagnetic  coupling  among  two-dimensional  electron  gas, 
graphene  micro-ribbon  and  surface  plasmons”,  a  self-consistent  theory  involving  Maxwell 
equations  and  a  density-matrix  linear-response  theory  is  solved.  The  conditions  of  the  solution 
include  an  electromagnetically-coupled  doped  graphene  micro-ribbon  array  and  a  quantum-well 
electron  gas  sitting  at  an  interface  between  a  half-space  of  air  and  another  half-space  of  a  doped 
semiconductor  substrate  which  supports  a  surface-plasmon  mode  in  our  system.  The  coupling 
between  a  spatially-modulated  total  electromagnetic  field  and  the  electron  dynamics  in  a  Dirac- 
cone  of  a  graphene  ribbon,  as  well  as  the  coupling  of  the  far-field  specular  and  near-field  higher- 
order  diffraction  modes,  is  included  in  the  derived  electron  optical-response  function.  Full 
analytical  expressions  are  obtained  with  non-locality  for  the  optical-response  functions  of  a  two- 
dimensional  electron  gas  and  a  graphene  layer  with  an  induced  bandgap,  and  are  employed  in  our 
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numerical  calculations  beyond  the  long-wavelength  limit.  Both  the  transmissivity  and  reflectivity 
spectra,  as  well  as  their  dependence  on  different  configurations  of  our  system  and  on  the  array 
period,  ribbon  width,  graphene  chemical  potential,  drift  velocity  of  quantum-well  electron  gas 
and  bandgap  in  graphene,  are  studied.  Moreover,  the  transmitted  E-field  intensity  distribution  is 
calculated  to  demonstrate  its  connection  to  the  mixing  of  specular  and  diffraction  modes  of  the 
total  electromagnetic  field.  An  externally-tunable  electromagnetic  coupling  among  the  surface, 
conventional  electron-gas  and  massless  graphene  intraband  plasmon  excitations  is  discovered 
and  explained.  Furthermore,  a  comparison  is  made  between  the  dependence  of  the  graphene  - 
plasmon  energy  on  the  ribbon  width  and  chemical  potential  in  this  paper  and  the  recent 
experimental  observation  given  by  Ju,  et  ah,  [Nature  Nanotechnology  6,  630  (2011)]  for  a 
graphene  micro-ribbon  array  in  the  terahertz-frequency  range. 

We  investigated  the  energy  transfer  properties  of  a  novel  nanocomposite  consisting  of  a  quantum 
dot  and  graphene  nanoparticle.  The  nanocomposite  is  embedded  in  a  photonic  crystal.  Energy 
transfer  occurs  when  electronic  excitons  in  the  quantum  dot  and  surface  plasmons  in  graphene 
interact  through  the  dipole-dipole  interaction.  Two  dressed  correlated  excitons  are  created  in  the 
quantum  dot,  corresponding  to  two  peaks  in  the  loss  function.  It  was  shown  that  the  exciton 
energies  may  be  transferred  to  graphene  through  the  dipole-dipole  interaction.  Additionally,  this 
energy  transfer  may  be  switched  on  and  off  by  applying  a  pump  laser  to  the  photonic  crystal  or 
by  adjusting  the  strength  of  the  dipole-dipole  interaction.  The  intensity  of  the  peaks  for  the  loss 
function  increases  when  the  number  of  graphene  layers  is  increased  or  the  separation  between  the 
quantum  dot  and  graphene  is  decreased.  The  principle  of  our  nanocomposite  may  be  employed  to 
fabricate  nano-biosensors,  optical  nano-switches,  energy  transfer  devices  and  quantum  tele¬ 
transportation  devices. 

4  Results  and  Discussion 

The  properties  of  the  high  energy  optical  ^--plasmons  of  simple  hexagonal  intrinsic  graphite 
were  calculated  within  the  self-consistent-field  approximation.  The  plasmon  frequeqpy  co 
was  detennined  as  functions  of  the  wave  vector  q\\  along  the  hexagonal  plane  in  the  Brillouin 
zone  and  the  perpendicular  component  q~.  When  qz  =  0,  the  plasmons  exist  for  wave  vectors 
larger  than  a  critical  wave  vector  qc  ~  0.02  A1 .  These  plasmons  are  isotropic  within  the  plane 
in  the  long  wavelength  limit  and  for  q\\  >  qc.  As  the  wave  vector  qy  increases,  the 
plasmon  frequency  strongly  depends  on  the  magnitude  and  direction  {(f>)  of  q  .  With 
increasing  <j),  the  dispersion  relation  of  co-(q\\)  is  gradually  changed  from  quadratic  to  nearly 
linear  fonn.  There  are  many  significant  differences  for  the  ^-plasmon  dispersion  relations 
between  2D  graphene  and  3D  dd -stacked  graphite.  They  include  threshold  frequency,  critical 
transferred  momentum,  q\\-  and  ^-dependence  and  ^--plasmon  bandwidth.  This  result  reveals 
that  interlayer  interactions  could  enhance  anisotropy  of  in-plane  ^--plasmons.  For  choseny  > 
qc,  we  also  obtain  the  co  as  a  function  of  qz  and  show  that  there  is  an  upper  bound  on  q-  for 
plasmons  to  exist  in  graphite. 

We  presented  the  formalism  and  numerical  results  for  the  energy  loss  of  a  charged  particle 
scattered  at  an  arbitrary  angle  from  epitaxially-grown  multilayer  graphene  (MLG).  It  was 
compared  with  that  of  free-standing  graphene  layers.  Specifically,  we  investigated  the  effect  of 
the  substrate-induced  energy  gap  on  one  of  the  layers.  The  gap  yields  collective  plasma 
oscillations  whose  characteristics  were  qualitatively  and  quantitatively  different  from  those 
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produced  by  Dirac  fermions  in  gapless  graphene.  The  range  of  wave  numbers  for  undamped  self- 
sustaining  plasmons  was  increased  as  the  gap  was  increased,  thereby  increasing  and  red-shifting 
the  MLG  stopping  power  for  some  range  of  charged  particle  velocity.  We  also  applied  our 
formalism  to  interpret  several  distinct  features  of  experimentally-obtained  electron  energy  loss 
spectroscopy  (EELS)  data. 

5  Conclusions 

We  calculated  the  edge-plasmon  excitation  spectrum  within  the  anti-crossing  bulk  bandgap 
region  for  an  inverted  HgTe/CdTe  quantum  well  by  employing  the  Bernevig,  Hugues  and  Zhang 
(BHZ)  model  within  the  random-phase  approximation  (RPA).  Our  proposed  model  system 
consists  of  a  single  component  electron  helical  liquid  (HL)  in  a  semi-infinite  quantum  well,  in 
which  a  topological  state  that  is  localized  around  the  system  edge  can  exist.  With  linearly- 
polarized  incident  light  as  a  perturbation  to  such  a  single  component  HL,  the  unique  charge 
dynamics  for  collective  excitation  of  these  edge-bound  electrons  with  broken  time-reversal 
symmetry  was  investigated.  The  plasmon  dispersion  cop(q)  of  the  single  component  HL  has  the 
form  cop(q)  ~  -a>oq  ln(qa)  in  the  long-wavelength  limit,  which  is  in  sharp  contrast  with  cop(q)  ~  - 
<x>oq  [-ln(qa)]  for  a  one-dimensional  electron  gas  (1DEG)  in  a  quantum-wire  system.  Here,  the 
plasmon  wave  number  in  a  conventional  1DEG  is  scaled  with  a  characteristic  width  W,  while 
that  in  our  model  system  is  scaled  with  a  lattice  constant  a.  Similar  to  interband  plasmons  in  a 
metallic  armchair  graphene  nanoribbon,  a>o  was  found  to  be  independent  of  the  linear  electron 
density  for  the  intraband  plasmon  in  our  system.  Besides  the  spin  factor  of  two  and  W  scaling  for 
the  wave  number,  the  dispersion  relation  for  the  collective  excitation  of  the  two-component  HL 
is  the  same  as  that  in  an  armchair  graphene  nanoribbon.  The  particle -hole  excitation  region  in  our 
system  is  found  to  be  collapsed  into  a  straight  line,  instead  of  a  wide  region  for  a  conventional 
1DEG.  The  plasmon  energies  of  the  two  and  single-component  HL  are  spectrally  separated 
although  the  same  particle-hole  excitation  region  is  shared  by  both  of  them. 

A  dilute  distribution  of  magnetic  impurities  was  assumed  to  be  present  in  doped  graphene.  We 
calculated  the  interaction  energy  between  two  magnetic  impurities  which  are  coupled  via  the 
indirect-exchange  or  Ruderman-Kittel-Kasuva-Yosida  (RKKY)  interaction  by  the  doped 
conduction  electrons.  Our  model  was  a  half-filled  AB-lattice  structure.  The  calculations  were 
based  on  the  retarded  lattice  Green's  function  formalism  in  momentum-energy  space  which  was 
employed  in  linear  response  theory  to  determine  the  magnetic  susceptibility  in  coordinate  space. 
Analytic  results  are  obtained  for  gapped  graphene  when  the  magnetic  impurities  were  placed  on 
the  A  and  B  sublattice  sites  of  the  structure.  This  interaction,  which  is  important  in  determining 
spin  ordering,  has  been  found  to  be  significantly  different  for  AA  and  BB  exchange  energies  in 
doped  graphene  due  to  the  existence  of  an  energy  gap,  and  is  attributed  to  a  consequence  of  the 
local  fields  not  being  equal  on  the  A  and  B  sublattices.  For  doped  graphene,  the  oscillations  of  all 
three  RKKY  interactions  from  ferromagnetic  to  antiferromagnetic  with  increasing  Fenni  energy 
is  significantly  modified  by  the  energy  gap  both  in  magnitude  and  phase.  Additionally,  the  AB 
exchange  energy  may  be  modified  by  the  presence  of  a  gap  for  undoped  graphene  but  not  for 
doped  graphene  due  to  the  dominance  of  doped  conduction  electrons. 
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